Trends in Elasticity and Electronic Structure of Transition-Metal Nitrides and 

Carbides from First Principles 
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The elastic properties of the Bi-structured transition-metal nitrides and their carbide counterparts 
are studied using the ab initio density functional perturbation theory. The linear response results of 
elastic constants are in excellent agreement with those obtained from numerical derivative methods, 
and are also consistent with measured data. We find the following trends; (1) Bulk moduli B and 
tetragonal shear moduli G' — (Cn — Ci2)/2, increase and lattice constants ao decrease rightward or 
downward on the Periodic Table for the metal component or if C is replaced by N; (2) The inequality 
B > G' > G > holds for G = C44; (3) G depends strongly on the number of valence electrons per 
unit cell (Zv)- From the fitted curve of G as a function of Zv, we can predict that MoN is unstable 
in Bi structure, and transition-metal carbonitrides {e.g. ZrCa:Ni_a:) and di-transition-metal carbides 
{e.g. Hfa:Tai_a:C) have maximum G at Zv ~ 8.3. 

PACS numbers: 62.20.Dc, 71. 15. Mb, 74.70.Ad 
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Elasticity describes the response of a crystal under ex- 
ternal strain and provides key information of the bonding 
strength between nearest-neighbor atoms. The informa- 
tion obtained from accurate calculation of elasticity is 
essential for understanding the macroscopic mechanical 
properties of solids and for the design of hard materials. 
Novifadays it is possible to calculate elasticity using ab 
initio quantum-mechanical techniques, and ab initio cal- 
culations have proven to be very powerful in not only pro- 
viding accurate elastic constants or moduli in good agree- 
ment with measurements but also predicting elastic- 
ity at extreme conditions of high temperatures and high 
pressures 0, ll{ , which are not easily accessible to exper- 
iment but have wide applications in the fields ranging 
from solid-state physics to seismology. Most previous ab 
initio calculations of elasticity used finite strain meth- 
ods within the framework of the density-functional the- 
ory (DFT). The development of density- functional per- 
turbation theory (DFPT) makes it possible now to obtain 
elastic constants directly and more accurately 0, Q . 

Transition-metal nitrides and carbides in the rocksalt 
{Bi) structure are widely used for cutting tools, mag- 
netic storage devices, generators and maglev trains due 
to their high hardness, high melting points and oxida- 
tion resistance 0. These excellent properties are as- 
sociated with their unusual electronic bonding. The 
relatively high superconducting transition temperature 
in some of these compounds, reaching nearly 18 K in 
NbCi_a:N2: j^, indicates a strong electron-phonon in- 
teraction. Many theoretical studies of their electronic 
structure 0, IE I1J> UM have revealed an unusual 

mixture of covalent, metallic, and ionic contributions to 
bonding which must ultimately lie at the root of their un- 
usual properties. Specifically, it was found that the 
hardness enhancement of these materials can be under- 
stood on a fundamental level in terms of their electronic 
band structures. But the general trends of elasticity and 



electronic structure among the transition-metal nitrides 
and carbides remain unclear and a challenge for engineer- 
ing hard materials. 

In this Letter we perform systematic calculations of 
elasticity of the selected transition-metal nitrides and 
carbides within the framework of the density functional 
perturbation theory. The accuracy of this method is also 
examined by numerical finite strain methods. We obtain 
interesting trends of the variation of the bulk and shear 
moduli among various compounds. We demonstrate the 
electronic origin of such trends in these materials. 

Elastic constants are defined as 
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— SijOki — Cijki—SijCfki , (1) 



where a is stress, e is strain, Q and Ed are the unit-cell 
volume and energy, S is the Kronecker delta function, and 
the Latin indices run from 1 to 3. At zero hydrostatic 
pressure Cijki = Cijki. One can apply a small strain and 
calculate the change of energy or stress to obtain elastic 
constants. Direct calculations of stress are possible from 
the quantum mechanical theory of stress |l4| . 

Alternatively elastic constants can be viewed as the 
linear response (LR) of an undisturbed system by homo- 
geneous strains (macroscopic distortions of the crystal), 
and can be directly computed by the DFPT How- 
ever, strain perturbation is much more difficult to cope 
with within the DFPT than atomic displacement pertur- 
bation (phonon) in that a homogeneous strain changes 
the boundary conditions of the Hamiltonian of an infi- 
nite system, and the DFPT requires the same basis set 
for the undisturbed and disturbed systems. 

Recently Hamann et al. ^ applied the reduced- 
coordinate metric tensor method to the LR of strain type 
perturbations. Their approach is based on an energy 
expression in terms of reduced coordinates, which are 
defined in both real and reciprocal space using primitive 
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FIG. 1: Theoretical (a) lattice constants ao, (b) bulk moduli 
B, and (c) zone-center band gap AEr for selected transition- 
metal nitrides and carbides with Bi structure. Filled symbols 
denote nitrides, and open symbols refer to carbides. 



direct (R) and reciprocal (G) lattice vectors. If the re- 
duced coordinates are denoted with a tilde, a real space 
vector X and a reciprocal space vector K are given by 



— ^ RijXj, Ki — ^ GijGj , 



(2) 



respectively. Essentially every energy term of the DFT 
can be expressed as dot products of vectors in real or 
reciprocal space, and 

X' . X ^ ^ X^S,,X, , K' • K = ^ k:T,,K, , (3) 

ij ij 

where the metric tensors a and T are 



(4) 



The DFT energy derivatives with respect to strain act 
only on the metric tensors, and the first order derivatives 
are 
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-GkiGij—GuGkj- (5) 



and second order derivatives can be directly obtained. 
After this transformation, strain has an equal footing 
with other coordinates, and the wave functions have in- 
variant boundary conditions, so that the strain deriva- 
tives such as elastic constants can be evaluated in a way 
similar to other derivatives. Note that the DFPT ex- 
pressions give Gijki, not Cijki- Elastic constants were 
computed for an insulator AlP [^, and the agreement 
between DFPT and numerical finite strain methods is 
perfect. Although the DFPT is directly applicable to 
metals, it needs special care for the Fermi surface due to 
partially occupied states. Here we examine the accuracy 
of DFPT for transition- metal nitrides and carbides. 

Our first principles calculations are based on the DFT 
within the local density approximation (LDA). We used 
ABINIT (version 4.3.3) 15], which uses plane-wave basis 
sets and norm-conserving pseudopotentials. The pseu- 
dopotentials of C, N, Zr, Nb, Hf and Ta were generated 



TABLE I: Calculated bulk moduli and elastic constants (GPa) 
of NbN in Bi structure. Three different approaches were used. 
See text for the meaning of theory 1, 2, and 3. 



method 


B 


Cii 


Cl2 


C44 


theory 1 


353.6 


742.6 


159.1 


76.4 


theory 2 


353.6 


739.4 


160.8 


75.5 


theory 3 


353.7 


738.9 


161.1 


74.8 



using the OPIUM program 16]. The cutoff energy of 
plane-wave basis is 50 Ha. A dense 16 x 16 x 16 k-point 
mesh was used over the Brillouin zone (BZ), and the cold 
smearing method [T^ was performed for BZ integration. 

Equilibrium lattice constants oq and bulk moduli B 
were evaluated from fitting energy and volume data to 
the Vinet equation of state (l^. Compared with the 
available experimental data [a, llSl l2Ji i "^o are about 
1.0-1.5% smaller, while in average B are 15% larger. 
These are typical LDA errors. The presence of N or C 
in transition-metal compounds leads to the hybridized p- 
d bonding, and its strength determines lattice constants 
and bulk moduli. As illustrated in Fig. QJa) and^Jb), 
ao decrease and B increase as the transition-metal pro- 
ceeds rightward (e.g. Zr Nb) or downward (e.g. Zr — > 
Hf) on the Periodic Table. Smaller ao and larger B in- 
dicate stronger bonding. The former trend is because of 
more filling of the metal d-bands, which enhances the d-p 
bonding 10]. The later trend is a result of the presence 
of /-electrons in core states, which repel d-orbitals out 
of core regions and in turn make them bond with C or 
N p-electrons more tightly. Available experimental data 
support these theoretical conclusions. Fig. also show 
that nitrides have smaller ao and larger B than their 
corresponding carbides, because N has one more valence 
electron than does C so that the d-p bonding strength in 
nitrides is higher than that in their carbide counterparts. 
Experiment and theory agree well with the trend of ao, 
but the existing measured moduli of NbN and ZrN are 
lower than those of NbC and ZrC, respectively. Because 
different experimental methods have different level of ac- 
curacy, and experimental circumstances {e.g. tempera- 
ture, purity, etc) could affect measured results greatly, 
experimentally re-scrutinizing B is needed to clarify the 
contradiction with theory. 

Electronic band structures were calculated at zero 
pressure, and it can be seen in Fig. El that all these 
nitrides and carbides in Bi structure are metals, since 
there is a direct band overlap at the X point. At the 
zone-center (F-point), there is a finite band gap AEy, 
except for NbC and TaC. AEr (Fig. IJc)), the d-band 
splitting breadth, increases as the transition-metal goes 
leftward on the Periodic Table due to an upward shift of 
the metal d-bands as a result of less filling ^3]j or goes 
downward because the presence of /-electrons promotes 
metal d-band energies. Carbides have smaller AEr than 



3 



their nitride counterparts since the carbon p-bands have 
higher energy than nitrogen. 

Elastic constants were computed in the Voigt notation, 
in which the Greek subscripts of C^^ run from 1 to 6. 
Crystals in cubic structure have only three independent 
elastic constants, namely Cn, C12, and C44, and the bulk 
modulus is 

B=i(Cn+2Ci2) . (6) 
One can apply the following strains 
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to distort the lattice vectors, R' = (1 + e)R. The result- 
ing changes of energy density {U — Ed/il) are associated 
with elastic constants, 

f/tctr = ^(Cii-Ci2)^'+0(^') , (7) 

C/orth = + 0(,54) , (8) 

respectively. Cn and C12 can be obtained from Eqs. 
© and 10, and C44 from Eq. ((HJ. We denote this fi- 
nite strain approach based on total energy as "theory 
1". Another way is to compute the stress tensor ele- 
ments for selected strains. For Cn and C12 we applied a 
strain tensor with en = S (zero for other elements), and 
Cii = Uii/S, C12 = a'22/S- For C44 the applied strain 
was ei2 = £21 — <5, and C44 ~ ai2/2S. This is called 
"theory 2". We used S = ±0.002, ±0.004. The third 
method (theory 3) is the DFPT of LR of strain pertur- 
bation For comparison, elastic constants of NbN 
obtained from these three methods are shown in Tabled 
where B were calculated from Cn and C12 by Eq. © 
for theory 2 and 3. High level of agreement between the 
DFPT and the energy (theory 1) or stress (theory2) was 
achieved, and similar agreement was also made for other 
compounds. It confirms the reliability of the LR theory 
of strain perturbation for elastic constants of metals. 

As summarized in Table Ull present DFPT elastic con- 
stants are in accordance with experimental results. We 
point out that the uncertainty in neutron scattering mea- 
surements of elastic constants could be as high as 10-15 
%, and the experimental data in Table Ull were measured 
at room temperature. The elastic stability criteria for a 
cubic crystal 22] at ambient conditions are 

Cn + 2Ci2 > 0, C44 > , and Cn - C12 > , (9) 

i.e. all the bulk B, shear C = C44, and tetragonal shear 
C = (Cn — Ci2)/2 moduli are positive. Our results 
satisfy all three criteria in Eq. @, and it follows that 
these materials in rocksalt structure are stable, consis- 
tent with experiment. In general _B > G > C, but for 
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FIG. 2: Electronic band structures of selected transition- 
metal nitrides and carbides with Bi structure. The red curves 
are band structures of NbN and NbC under a finite shear 
strain (ei2 = £21 = 0.1). 



these materials both theory and experiment hold that 
B > C > G > 0, so the shear modulus G is the main 
constraint on stability. 

The tetragonal shear modulus G' measures the re- 
sponse of a crystal under volume conserving tetragonal 
shear strain, which involves stretching of metal-N or -C 
bonds and bending of metal-metal bonds. Because in 
these crystals Cn ^ C12, and Cn is determined by the 
nearest- neighbor interaction, similar to bulk modulus B, 
C have the same trends as B, as seen in Fig. Ela) that 
(1) G' increases as the metal goes rightward or down- 
ward on the Periodic Table; (2) Nitrides have bigger G' 
than their carbide counterparts. These trends generally 
are supported by experiment, and they are a result of 
the enhancement of d-p bonding due to adding valence 
electrons or the presence of /-electrons, as discussed in 
previous paragraphs for B. However, farther right on 
the Periodic Table {e.g. Nb Mo) valence electrons sat- 
urate the bonding states and begin to fill the antibonding 
states, so B and G' will go down after reaching a peak. 

The shear modulus G (C44) is the most important pa- 
rameter governing indentation hardness. The hardness 
of a material is defined as its resistance to another mate- 
rial penetrating its surface, and it is determined by the 
mobility of dislocations. In covalent hard materials, the 
bond-breaking energy under plastic deformation and the 
bond-restoring energy under elastic shear strain are very 
similar, so that one of the determining factors of hard- 
ness is the response of covalent bonds to shear strain ■ 
We classified these materials according to their number 
of valence electrons per unit cell {Zy) in Fig. ^h), which 
shows that G is intimately related to Zy ■ At Zy = 8 and 
9, these crystals have large G, but a modest decrease oc- 
curs when Zv rises from 8 to 9. Further increasing Zy 
to 10 significantly lowers G. To explain this fascinating 
phenomenon, we investigated electronic band structures 
under finite shear strain. As an example, in Fig. |21the 
band structures (red dotted curves) of NbN and ZrC un- 
der a strain of ei2 = £21 =0.1 are displayed against zero- 
strain curves. Under shear strain, p-d hybridized orbitals 
split at F-point. The energy of the fourth valence band 
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FIG. 3: Theoretical (a) tetragonal shear moduh G' — (Cn — 
Ci2)/2, and (b) shear moduli G = C44 with respect to number 
of valence electrons per unit cell (Zv)- The solid curve in (b) 
is obtained by quadratic fitting to average shear moduli at 
each Zv- 



TABLE II: Elastic constants (GPa) of selected transition- 
metal nitrides and carbides with Bi structure. Here theoret- 
ical results were obtained from the DFPT, and experimental 
data are in parenthesis. 



system 


Cii 


C12 


C44 


ZrN 


616 (471") 
522 (470') 
739 (608") 
667 (620*") 
694 (679") 
574 (500'') 


117 (88") 


130 (138") 


ZrC 
NbN 


110 (100'') 

161 (134") 


160 (160'') 
75 (117") 


NbC 
HfN 


163 (200'') 
112 (119") 


161 (150'') 
135 (150") 


HfC 


107 (114") 


180 (180'') 


TaN 


783 


167 


20 


TaC 


740 (550'') 


165 (150'') 


176 (190'') 



"Neutron scattering measurements, Ref. llflL 
''Estimation from phonon dispe rsio n curves, Ref. |22 
'^Ultrasonic measurements, Ref. l2ll 



increases dramatically (red thick lines in Fig. [SJl in the 
L-F-K section, and the energy of the fifth band decreases. 
Strain £4 involves shearing metal-N or -C bonds, and the 
direction-sensitive bonding character results in opposite 
movements of these two bands. For Zy — 8 (see band 
structures of ZrC and HfC in Fig. along the L-F-K 
hne, the fourth valence band are nearly fully occupied 
while the fifth band is empty, so a large G is expected. 
On the other hand, at Zy = 9 (NbC, TaC, ZrN, and 
HfN), valence electrons begin to fill the fifth band, which 
leads to a negative contribution to G. However, because 
the occupation on the fifth band for NbC and TaC is 
tiny, G of NbC and TaC are very close to those of ZrC 
and HfC, respectively. But ZrN and HfN have more elec- 
trons filling the fifth valence band around F-point than 
do NbC and TaC, so G of ZrN and HfN are noticeably 
smaller than those of ZrC and HfC, respectively. A steep 
decline of G occurs for Zy = 10 (NbN and TaN) because 
of substantial filling on the fifth valence band. Note that 
TaN has the largest B and G' among these materials, but 
its G is the smallest; and interestingly, TaC have large 
values of all B, G' and G. 

We fitted three average G points to a simple quadratic 
form and extrapolated the fitted curve to Zy — 7 and 
11 (Fig. Elb)). This fitted curve predicts a G w 105 
GPa at Zy = 7, and very close to the actual calculated 
G = 114 GPa for YC. YC has smaller G than ZrC due 
to its half filling of the fourth valence band. Our the- 
ory predicts a negative G at Zy = 11, and indeed we 
computed G44 = -67 GPa for MoN. Thus MoN is un- 
stable in rocksalt structure because extra electrons fill 
to the antibonding states, and the energy gain from the 
fourth band is less than the energy loss from the fifth 
band under shear strain. Previous first-principles cal- 
culations drew the same conclusion. We note that 
the peak of the fitted curve in Fig. |3fb) is roughly at 
Zy = 8.3, which corresponds to full occupation of the 
fourth valence band without filling on the fifth band, and 
an optimum G is expected. Non-integer Zy can be real- 



ized by transition-metal carbonitrides such as ZrCxNi-^ 
or di-transition-metal carbides such as ZrajNbi-^C Our 
estimation of Zy w 8.3 for maximum G agrees surpris- 
ingly well with experiment 0, 123I [2^ of Zy = 8.4 and 
the ab initio virtual crystal method 111 of Zy in a 
range of 8.3-8.5. Since both TaC and HfC have large 
G, we propose that Hfj.Tai_a,C with x « 0.3 has the 
largest shear modulus among these transition-metal car- 
bonitrides or di-transition-metal carbides in Bi struc- 
ture. Zra;Nbi_a;C and Tia;Vi_a;C should have compa- 
rable G values to Hfa;Tai_2;C, but they have advantages 
of much less weight. 

We have carried out first-principles calculations of 
structural and elastic properties of selected transition- 
metal nitrides and carbides. We computed elastic con- 
stants of metals using the density-functional perturba- 
tion theory, and our results proved the accuracy of the 
metric tensor approach Q of the linear response theory 
of strain perturbation. It has been demonstrated that 
these crystals have similar trends for bulk B and tetrag- 
onal shear G' moduli, but different trend for shear mod- 
uli G. B and G' are associated with strains essentially 
keeping the crystal symmetry, and all valence bands re- 
spond uniformly by moving up or moving down. Shear 
strain £4 lowers the crystal symmetry significantly, and 
it causes p-d orbitals splitting at zone-center and it has 
opposite effects on the fourth and fifth valence bands. 
Such a systematic study can help clarify the ambiguity 
rising from different experimental methods, and predict 
new materials with better properties. 

This work was supported by the Office of Naval 
Research (ONR) Grant number N00014-02-1-0506, the 
Department of Energy (DOE) Grant Nos. DEFG02- 
02ER4595 and DEFC03-03NA00144, and DOE ASC 
subcontract B341492 to Caltech DOE W-7405-ENG- 
48(REC). ZGW thanks the helpful discussions with E. 
J. Walter. 



5 



[1] C. S. Zha, H. K. Mao, and R. J. Hemley, Proc. Natl. 
Acad. Sci. U.S.A. 97, 13494 (1999). 

[2] L. Stixrude and R. E. Cohen, Science 267, 1972 (1995). 

[3] G. Steinle-Neumann, L. Stixrude, and R. E. Cohen, Na- 
ture (London) 413, 57 (2001). 

[4] S. Baroni, S. de Gironcoli, A. Dal Corso, and R. Gian- 
nozzi. Rev. Mod. Phys. 73, 515 (2001). 

[5] D. R. Hamann, X. Wu, K. M. Rabe, and D. Vanderbilt, 
arXiv:cond-mat/0409269 

[6] L. E. Toth, Transition metal carbides and nitrides (Aca- 
demic Press, New York, 1971). 

[7] D. A. Papaconstantopoulos, W. E. Pickett, B. M. Klein, 
and L. L. Boyer, Phys. Rev. B 31, 752 (1985). 

[8] B. T. Matthias, Phys. Rev. 92, 874 (1953). 

[9] J. Chen, L. L. Boyer, H. Krakauer, and M. J. Mehl, Phys. 
Rev. B 37, 3295 (1988). 
[10] J. C. Grossman, A. Mizel, M. Cote, M. L. Cohen, and S. 

G. Louie, Phys. Rev. B 60, 6343 (1999). 
[11] S.-H. Jhi, J. Ihm, S. G. Louie, and M. L. Cohen, Nature 

(London) 399, 132 (1999). 
[12] C. Stampfl, W. Mannstadt, R. Asahi, and A. J. Freeman, 

Phys. Rev. B 63, 155106 (2001). 
[13] S.-H. Jhi, S. G. Louie, M. L. Cohen, and J. Ihm, Phys. 
Rev. Lett. 86, 3348 (2001). 



[14] O. H. Nielsen and R. M. Martin, Phys. Rev. Lett. 50, 
697 (1983); Phys. Rev. B 32, 3780 (1985). 

[15] X. Gonze, J.-M. Beuken, R. Caracas, F. Detraux, M. 
Fuchs, G.-M. Rignanese, L. Sindic, M. Verstraete, G. 
Zerah, F. JoUet, M. Torrent, A. Roy, M. Mikami, Ph. 
Ghosez, J.-Y. Raty, and D. C. Allan, Comput. Mater. 
Sci. 25, 478 (2002). 

[16] A. M. Rappe, K. M. Rabe, E. Kaxiras, and J. D. 
Joannopoulos, Phys. Rev. B 41, 1227 (1990); L Grin- 
berg, N. J. Ramer, and A. M. Rappe, Phys. Rev. B 62, 
2311 (2000). 

[17] N. Marzari, Ph.D. thesis, Cambridge University, 1996. 
[18] R. E. Cohen, C. Gulseren, and R. J. Hemley, American 

Mineralogist 85, 338 (2000). 
[19] X.-J. Chen, V. V. Struzhkin, Z. Wu, M. Somayazulu, J. 

Qian, S. Kung, A. N. Christens, Y. Zhao, R. E. Cohen, 

H.-K. Mao and R. J. Hemley, unpublished. 
[20] W. Weber, Phys. Rev. B 8, 5082 (1973). 
[21] H. L. Brown, P. E. Armstrong, and C. P. Kempter, J. 

Chem. Phys. 45, 547 (1966). 
[22] J. Wang, S. Yip, S. R. PhiUpot, and D. Wolf, Phys. Rev. 

Lett. 71, 4182 (1993). 
[23] M. HoUeck, J. Vac. Sci. Technol. A 4, 2661 (1986). 
[24] V. Richter, A. Beger, J. Drobniewski, L Endler, and E. 

Wolf, Mater. Sci. Eng. A 209, 353 (1996). 



